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The eigenvalues of the Hyperspherical billiard are calculated in the semiclassical approximation. 
The eigenvalues where this approximation fails are identified and found to be related to caustics 
that approach the wall of the billiard. The fraction of energy levels for which the semiclassical 
error is larger than some given value is calculated analytically (and tested numerically) and found 
to be independent of energy. The implications for other systems, in particular integrable ones, are 
0^ ' discussed. 
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I. INTRODUCTION 



Exact solutions to physical problems are rare. In most cases one has to resort to approximate 
solutions. In quantum mechanics the Wentzel-Kramers-Brillouin (WKB) method, along with per- 
turbation theory, is probably the most common method used to obtain approximate solutions, (for 
reviews see 0^0] )• The WKB method is a formal h expansion for the wave function, that expresses 
its rapid oscillations in the semiclassical limit. Using this expansion combined with an appropriate 
boundary condition results in an approximate quantization condition. Therefore, high order approx- 
q • imations for the wave function can be used to improve the accuracy of the eigenvalues. A systematic 
way to obtain approximate quantum eigenvalues using WKB was developed first by Dunham [|J, 
was improved by Bender and coworkers || and is summarized in 0. In some cases the resulting 
series for the eigenvalues converge to the exact ones, but generally the resulting expansion will be 
^ ■ an asymptotic series. 

A surprisingly small number of papers were devoted to the accuracy of semiclassical methods, 
c$ ', and even less for high orders or high dimensional systems. A naive estimate for the accuracy of 
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semiclassical methods is that the leading semiclassical approximation is accurate to order h and 



therefore the resulting error is of order h 2 . For example, substitution of the Van Vleck propagator 
into Schrodingers equation will not solve the equation exactly, and a reminder of order h 2 will 
result @. The mean level spacing scales like h D , where D is the dimension of the system, thus the 
relative error scales like h D ~ 2 and the semiclassical method fails for high dimensional systems. The 
prefactor of the semiclassical error is important (in particular for D = 2), since individual energies 
can be found semiclassically only if the error is less than the mean level spacing. For integrable 
systems the WKB expansion enables one to find analytic estimates for the error in the energies 
and this will be demonstrated for a class of such systems in this paper. In addition it may shed 
light also on the accuracy of the semiclassical approximation for other dynamical systems. The 
common argument for the failure of the semiclassical approximation for chaotic systems in D > 2 
and possibly for D = 2 was recently challenged by Primack and Smilansky ||. For chaotic systems 
energy is the only quantum number and this is the only identity of a level. The mean level spacing 
is determined by the smooth Weyl term. This term can be determined for billiards to a high order 
in Plancks constant H. Primack and Smilansky introduced a method to evaluate the error in the 
energy levels resulting from the fact that oscillatory terms are known only to the leading order in H. 
They concluded that the error in the evaluation of single levels compared to the mean level spacing 
diverges at most as | In h\ in the limit h — > 0. Dahlqvist B estimated effects of diffraction for the 
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D = 2 Sinai billiard and concluded that ignoring diffraction, as is done in the standard semiclassical 
approximation, may result in errors at least of the order of the mean level spacing. 

Here we will study the much simpler case of integrable systems when an expression for the energy 
levels can be obtained to arbitrary order in h and the levels that are poorly approximated in the 
leading order of the semiclassical approximation can be characterized. An estimate for the error 
of the semiclassical quantization for the Hyperspherical billiard, that is a generalization of the 
circle billiard, will be obtained. This billiard has been subjected to several works regarding the 
semiclassical accuracy. Prosen and Robnik explored the error for the energies of the circle billiard 
numerically [TO] . They found that the mean error increases with the energy, and concluded that 



the semiclassical approximation fails in this simple system. Boasman used the boundary integral 
method to find the quantization error for several billiards, including the circle billiard ||11|| . The 
error was found to be a small fraction of the mean level spacing for general two dimensional billiards. 
For the circle billiard he found that for most of the energies the semiclassical error is a constant, 
which is a small fraction of the mean level spacing. There were large errors for some eigenenergies, 
the corresponding eigenfunctions were found to be affected by caustics, and it was stated that the 
fraction of these poorly approximated states decreases with energy. 

Following earlier work of Agam |12| we use a certain WKB expansion, keeping the classical quanti- 
ties fixed, that gives identical quantization condition to the one obtained from the Debye expansion 
of the Bessel function. This expansion enables us to estimate the fraction of states where the error 
in the leading semiclassical approximation exceeds some value. The relation to caustics turns out 
transparent. For integrable systems this is a better measure then the mean error since some eigenen- 
ergies, that can be clearly characterized, may have extremely large errors. The reason is that for 
integrable systems, contrary to chaotic ones, there are other quantum numbers in addition to the 
energy. These help to characterize the groups of quantum eigenstates that are badly approximated 
in the leading order of the semiclassical approximation. 

The semiclassical expansion for the eigenvalues of the Hyperspherical billiard to second order in 
h is derived in Sec. 01 with the help of the standard semiclassical expansion, as well as the Debye 
expansion of the Bessel functions. In Sec. [ITJ the density of eigenvalues where the semiclassical 
approximation is poor (a term that is defined in that section) is derived analytically and tested 
numerically. The results are summarized in Sec. [TV] and the implications for other systems are 
discussed. 



II. THE WKB EXPANSION FOR THE HYPERSPHERICAL BILLIARD 



In this section the semiclassical quantization for the Hyperspherical billiard will be developed up 
to second order in h. First a short introduction to a systematic WKB expansion for one dimensional 
problems will be given. Then the Debye expansion for the Bessel function will be used to derive the 
second order quantization condition. This derivation is much simpler than the derivation that uses 
a systematic WKB expansion. In order to verify that the quantization condition derived from the 
Debye expansion of the Bessel function is indeed a result of a semiclassical expansion a WKB series 
will be developed for the radial wave function. Quantization of this series indeed gives the same 
quantization condition that was obtained using the Debye series (at least up to second order). 

A scheme for quantization using high order WKB expansion is known, at least for an ordinary 
differential equation such as the equations found for separable Hamiltonians 0-^,13]. For one di- 
mensional systems the Schrodinger equation is 



2m dx 2 



+ V(x) 



ip{x) = Eip(x) 
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The wave function is written in the form 



il>(x) = A exp \j:<r(x)) , (2) 
and the phase is expanded in powers of % 

while A is a normalization constant. The Schrodinger equation can be solved order by order in h 
leading to a recursion relation for the expansion functions [f§-|6|] 

a'l = 2m(E - V(x)) 

n 

E «-k + <-i = o- (4) 

fc=0 

The quantization condition is derived by demanding that the wave function is single valued 

I — ] j> dak = 2nhn. (5) 
fc=0 V 1 J J 

The first odd term <j\ has the form of a logarithmic derivative, do\ = — jd^na'^). Each turning 
point gives a simple zero of <7q 2 , and the contour integral for <Ji counts these zeros and results in the 
Maslov index. All of the other odd terms are real, and usually without cuts in the complex x plane. 
Thus these terms will not contribute to the quantization condition [||]. Therefore, the quantization 
condition becomes 




j> da 2k = 27fh(n + m/4) (6) 



where m is the Maslov index to be evaluated later (for the simple problem with two turning points, 
m is just 2, the number of zeros of <r 2 ) . The fact that only even terms contribute to the quantization 
makes this quantization method efficient since the correction to the eigenvalues will be smaller by 
h 2 instead of h. This will motivate us transform the radial differential equation to a form without a 
first derivative. 

The specific problem that will be studied in this work is the D dimensional Hyperspherical billiard, 
namely a free particle inside a D dimensional ball (R 2 > Y^iLi X D with Dirichlet boundary conditions. 
Generalized spherical coordinates will be most convenient to solve the Schrodinger equation. The D 
dimensional Hamiltonian reduces to the Laplacian operator 

A \ 

(7) 

with the boundary condition ip{R) = on the wave functions. The generalized angular momentum 
operator is A s and its eigenvalues are — 1(1 + D — 2) [Q. The radial and angular variables can be 




separated as \l/ = i2(r)x(fi). For the Hyperspherical billiard the radial Schrodinger equation is 



Tl"{r) + —TZ\r) - -i ! -K{r) + ~^-^(r) = (8) 
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The exact solution can be easily obtained using the fact that the radial equation (|8|) can be reduced 
to the Bessel equation of order u, namely 



u 2 J" + uj' + (u 2 - v 2 )J = 



(9) 



where R{u) 



u 2 



' J{u), v = I H — 2~ an d u 



2mE 
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of the first kind of order z/, and the quantization condition is: 



\-r. The solutions are J v , the Bessel functions 



l2mE 



h 2 



R 



(10) 



This result is exact. The dimensionless variables 
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D 



ft 2 

will be introduced for simplicity It is obvious that z > 1. 

In order to make contact with the semiclassical expansion it is useful to introduce a large order 
expansion of the Bessel function since for small h, the order v is large. This is the Debye asymptotic 
expansion. 



KVy/ Z 2 — 1 



<x> 

cosC ■ 

m=0 



■l) m r(2m+|)q 2r , 



2m 



+ 



oc (_i)m r ( 2m+ |) a 
+ smC 2^ 



2m+l 



i, <:-n ^(2) 

— f and a = 1 , while % 
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J, and ct2 = Y28 + 
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Where ( = v z 2 — 1 — arccos 

385 1 irrai 
3456 ( 2 2-l) 2 ' " ' U±al- 

Semiclassical quantization, using the Debye expansion, is performed in orders of ft, by requiring 
that the wave function will have a zero at r = R, order by order in K Thus to the second order the 
quantization condition is: 



cosC + 



v^TIV8 24z 2 -l 



sin C = 0. 



(13) 



First we solve for the first order condition and correct it perturbatively. Define Co — C(zo) as the 
leading order solution. It satisfies 



cos Co = 0. 



(14) 



The second order solution can be expanded around the first order one: z = Zq + 5z. The term 
5( = ((zq + Sz) — C(zq) is of the order i and therefore small. Substituting in (|13|) one finds 



>y/z^T 



1 5 

77 + 



24 z 2 - 1 



(15) 



This last expression can be given in terms of the first order solutions Zq or the second order ones 
since the difference is of higher order in ~. The last step is to rewrite the first order condition as 
cos Co = cos(C — S() = 0. This results in 
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7r [n + 



(16) 



This is the second order semiclassical quantization condition. For every pair of quantum numbers 
(n, v) one can solve for z and obtain E sc (n,u). Since the quantization condition is determined by 
only two quantum numbers, the problem behaves like a system of two degrees of freedom , if one 
chooses to ignore degeneracies. The reason for denoting this quantization condition as a second 
order one will become clear shortly. 

We turn now to apply the systematic semiclassical expansion to (^j). Our goal is to show that the 
quantization condition that was obtained from asymptotic expansion of Bessel functions coincides 
with the one obtained from a WKB expansion. This equation has a first order derivative and 
therefore does not have the form of the one dimensional Schrodinger equation ([!]). To exploit the 
formal WKB expansion described above one has to eliminate the first order derivative from the 
equation. Substitution of 



K{r) 



in (H) leads to the following equation for 



— < 

2m 



[r + 



h 2 



2mr 2 



(I 



l-D 

r 2 i 



D-2, 



(17) 



<f>(r) = E(f){r). 



Thus after this substitution the resulting equation will be similar to that of a system with one 
degree of freedom and the WKB method will give better results. The angular momentum that 
appears in this equation is the semiclassical one (see (|TTD and App. [A]), as opposed to the exact 



angular momentum hJl(l + D — 2) as will be explained in what follows. These differ by an extra 
term that vanishes in the limit h — > 0. That gives one some freedom in defining the semiclassical 
variables. Thus, for example, we can define the limit as h — > while E and L 2 = h 2 l(l + D — 2) 
are constants. This prescription for the semiclassical limit is problematic, since the form of the 
resulting leading order wave function near zero and at infinity differs than from the one of the exact 
wave function. This problem is known for the spherical case (D = 3), where the exact angular 
momentum is replaced by the semiclassical one, a modification introduced by Langer [l§Jl|]- For 
the D dimensional generalized angular momentum the analog is the replacement of (/ +^y^) 2 — \ 
by (I + "^y 2 ") 2 - Therefore, we define the semiclassical limit as the limit where % — ► while E and 
L sc = h(l + ^y 2- ) are constant. The difference between these definitions of the angular momentum 



vanishes in the limit h 
equation 



0. As a result, we will apply the WKB approximation to the following 



h 2 (f)"(r) = 2m E 



T 2 h 2 ' 
^sc ~ 

2mr 2 



(19) 



Using the expansion (|3p for the wave function we find that the recursion relation (|j) changes to 
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(20) 



The first order contribution to the quantization is obtained from the O'th order of (|20] 
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<j' dr 
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2mEr 2 — L 2 c dr 
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(21) 



The variable z (see Eq. QTTJ) ) is the ratio between the radius of the billiard and the minimal distance 
of the classical orbit from the origin, or the radius of the caustic. Therefore, if z is close to 1 then 
the state is more affected by caustics and the semiclassical approximation deteriorates. The reason 
is that the wave function of these states is close both to the caustic and to the hard wall where the 
pure semiclassical approximation fails. The next term will give us the Maslov index. The integration 

over <j[ will give an extra phase of | due to the turning point at r = (resulting in a zero of 

Uq 2 as discussed after (||)). We have to include also the hard wall at r = R . To do this correctly one 
has to build the wavefunction for the sum of to solutions in such way that they cancel each other 
at the wall. The result is that the reflected wave acquires an extra phase of ir, therefore the total 
Maslov index is just 3. Consequently the quantization condition in the leading order is 
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(22) 



In App. [A] it is derived in the framework of the "old quantum theory" by direct action quantization. 

Our goal is now to obtain the second order semiclassical condition. A straight forward calculation 
leads to 
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2mEr 2 - L 2 SC 



(23) 



The contour integral over a' 2 is calculated in App. [FJ, and the result is: 
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Using the quantization condition (B|) up to the second order 
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and (|2l|) one obtains the second order quantization condition 
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(24) 



(25) 



(26) 



It is identical to (|16|), obtained from the Debye expansion. If L sc is kept fixed in the semiclassical 
limit then the limit 7i — ^ is equivalent to v —>■ oo. Note that if z is not too close to 1 the first term 
is much larger then the second one and dominates the result. 

It was demonstrated explicitly that the Debye series and the WKB method lead to the same 
quantization condition up to the second order. Furthermore, these expansions are valid in the same 
limit. Therefore we expect that these expansions give identical results in any order. This enables 
one to use the Debye asymptotic expansion, in this case, to obtain estimates for the semiclassical 
error. This will be the subject of the next section. 

Before we continue, it is worthwhile discuss briefly the semiclassical angular momentum L sc . In 
the App. ^ we calculated the semiclassical eigenvalues of the generalized angular momentum using 
first order WKB. One may ask whether it is possible to obtain the exact form of the angular 
momentum eigenvalues from a WKB expansion. It turns out that this is indeed the case. The WKB 
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approximation for the 3 dimensional angular momentum was treated by Robnik and Salasnich |fL7 
They computed the first several terms in the series for the angular momentum eigenvalues, and 
conjectured the general form of this series. This conjecture was later proved by Salasnich and 
Sattin ||18|| . The generalization for high dimensional systems is not difficult and requires only minor 
changes in the differential equation and the WKB series. The details are presented in App. 0. The 
convergent series leads to the correct eigenvalues for the generalized angular momentum (namely 

fi^Jl(l + D — 2)), while the first order in the WKB expansion is L sc — h(l + ^y^)- 

III. THE SEMICLASSICAL ERROR FOR HYPERSPHERICAL BILLIARDS 

The leading semiclassical eigenvalues correspond to the zeros of J^{yz\ the leading term in the 
Debye asymptotic expansion. Here J^iyz) denotes the i'th order term in this expansion. Let 

r(°)( 

■7755 ||lT| , [l2| , where j'J^iyz) is the first derivative of J^°\vz) with respect to z and from ( |i~2"D 

one finds, Az = jj^r where z is the leading order solution satisfying J'^(uzo) = (see ( |T4"[ ) 

\ dz J 

and (|T5|) . The result can be obtained also directly from these equations). 

AE 
2E 



Az = z(E ex ) — z(E sc ) } then to the leading order in this asymptotic expansion Az = — Jl ' " " ' 



If the difference Az is small then — = ^p. Using fllT]) 



\AE\ 



h 2 z 2 rl 5 1 



mR 2 z 2 — 1 



24 z 2 - 1 



(27) 



that is correct in the leading order in h. Note that the error AE is independent of the number 
of degrees of freedom D. In the limit z — > 1 the expression diverges. This divergence is related to 
fact that in this limit the caustic approaches the wall and the classical trajectories on the quantized 
torus are always adjacent to the caustic and to the hard wall. The approximation is best in the 
limit z — > 00. In this limit the angular momentum contribution to the energy is negligible, and the 
problem is effectively one dimensional. 

To estimate the accuracy of the approximation it may be more meaningful to measure the error in 
units of the mean level spacing A. In the D dimensional Hyperspherical billiard the leading order 
in the mean level spacing is given by the Weyl formula 

(2n) D h D 1 , . 

A ~ i - 2%) 

mR 2 V%DL°- 2 z D - 2 1 J 

D 

where Vb = /"""^p, is the volume of the D dimensional sphere of unit radius.. With ( p7|) and ( ^8|) 

1 (1+ 2 J 

the semiclassical error is: 



|A£| DL D ~ 2 z D 

1 1 be 



r 2 (i + f)2 s z 2 -i 



1 5 1 
^ + 



8 24 z 2 - 1 



h D ~ 2 



(29) 



In the calculation of the mean level spacing A, the levels were weighted with their degeneracy, 

\ AF ( 

leading to a large level density, resulting in a very small mean level spacing. The divergence of J -^- L 
for D > 2 in the limit h — > ( L sc and z fixed) and z — > 00 ( L sc and H fixed) is a result of the rapid 
increase of the number of levels in these limits. 

The Hyperspherical billiard exhibits a high degree of symmetry. Therefore there are only two 
quantum numbers, (n, I), that determine the quantization condition. As a result, if levels are not 
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weighted by their degeneracies the density of the degenerate energy levels is quasi two dimensional. 
Consequently, if the number of states that are affected by caustics is small than one can use WKB 
to determine single energy levels. Knowledge of the relative number of states that are badly approx- 
imated is crucial here. If such states are common than the semiclassical approximation will fail. In 
what follows the fraction of states that are badly approximated by WKB will be estimated. 

One can decompose the total spectrum using the angular momentum quantum number. Each 
sub-spectrum is a series of eigenvalues that depend on the quantum number n. The density of these 
states with respect to the energy can be approximated from the quantization condition fl22|) 



Pn - ^ 




V 



2txE\ e(y) 



E-e(u) , . 

(30) 



where e(u) = ^ R2 . The decomposed spectrum is not very interesting, since our purpose is to find 
what is the total fraction of states that are baldy approximated irrespective of the other quantum 
numbers. For this, the spectrum can be constructed by adding the contributions from different 
angular momenta. It is important to understand that when we sum over v the energy is kept fixed, 
thus we actually sum over the different quantum states (u, n) with energy near E. If the sum is over 
all the possible values of v then one is counting all the possible states and will get the total density 
of states. If one wants to sum only the states that are badly approximated one should sum only the 
states that satisfy the condition 1 < z < 1/a, where a is a measure for the maximally allowed error. 
From (|27|) it is clear that \AE\ is monotonically decreasing with z. To find the density of badly 
approximated states one sums only values of v satisfying this inequality. Replacing the restricted 
sum by an integral and changing variables to e, leads to 



E-e(u) mR 2 f E lE-ede 

(31) 



'2ttE^ e(u) 2nh 2 J^E ] i e E 
The resulting density of badly approximated states is 

mR 2 



P 



2nh 2 



arccos a — cr 



VT^x 2 ■ (32) 



The result depends only on a, which is a measure of the error but does not depend on energy. The 
density of badly approximated states does not change with the energy!! Thus, the probability that 
a state is badly approximated by the WKB method does not depend upon the energy of the state. 
It is important to notice that this expression is useful only in the limit a << 1. When a ~ 1 the 
connection between the allowed error and a is not accurate since the various terms in the Debye 
series are of comparable magnitude and the expansion cannot be terminated. Indeed, when all of the 
states are included, (for a = any state is considered to be badly approximated), the density of all 
the states is as expected since this is the total density of states according to the Weyl formula 
in two dimensions (the double degeneracy of states with positive and negative angular momentum 
quantum numbers is ignored here). 

To verify these results we computed the first 125336 exact and semiclassical levels in the D=4 
dimensional Hyperspherical billiard numerically. A convenient set of units is 2m — H — R — 1. The 
error in each eigenvalue was computed directly from comparison between the exact energy levels 
computed from (^) and their semiclassical approximation computed from (|2~2"D . A set of a values 
was chosen. Each value of a defines a semiclassical error by (|27D , with z = -. Then, the density of 
states with error larger than this semiclassical error was computed by counting the number of such 



8 



states in an energy interval which is much larger than the mean level spacing but small enough not to 
smear any possible energy dependence. This number is just the density of badly approximated states 
times the width of the energy interval. For each a the density of badly approximated states was 
found to be constant when the energy was changed, as expected from (0). The value of this density 
was compared to the density predicted by ([32]), and is presented in Fig. [I]. The approximation ([32]) 
was found to be excellent. 



0.20 



0.15 



a. 0.10 - 



0.05 - 



0.00 



0.0 




FIG. 1. The density of badly approximated states, (solid line - the approximation, circles - numerical results) 

It was found that using the fact that the system is integrable and the knowledge of its degeneracy, 
we can treat the spectrum of the Hyperspherical billiard as two dimensional. The density of badly 
approximated states was calculated and was found to be energy independent contrary to a previous 
statements ||11|| . Only in the limit where the angular momentum is negligible there are no large 
errors in quantization, since the system is effectively one dimensional. 



IV. DISCUSSION AND CONCLUSIONS 



The accuracy of the semiclassical approximation was studied for the Hyperspherical billiard. The 
approximation was found to fail for energies corresponding to states that are localized on such tori, 
that the caustic is close to the wall. As the caustic approaches the wall the region in space where 
the semiclassical approximation for the wave function is justified shrinks, and the quality of this 
approximation deteriorates. The semiclassical approximation for eigenvalues can be arbitrary bad 
(for z ~ 1, in our case). For this reason the average semiclassical error is meaningless since it may 
be dominated by few large contributions. It is much more reasonable to calculate the fraction of 
badly approximated states, namely the fraction of states where the semiclassical error exceeds some 
value. The density of such states was found to depend only on one parameter that measures the 
ratio between the angular momentum and its maximal possible classical value for a given energy 
(the parameter was denoted by a, and the regime where the semiclassical approximation fails is 
1 < z < -). The semiclassical error (]27|) depends only on z and not on the energy. Therefore we 
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concluded that the density of poorly approximated eigenvalues characterized by a is independent 
of energy. The semiclassical approximation does not improve as the energy increases, contrary to 



common belief (stated for example in Lazutkin found that convex billiards have tori of the 

KAM type along the boundary. He also proved that it is possible to use a WKB like method to 
obtain approximate solutions and energy levels on these tori JES) . For these states there is a caustic 
near a hard wall and we expect that the energy levels will be badly approximated. 

For integrable systems there are symmetries and the situation where degeneracies exist is quite 
typical. Since the symmetries are respected by the semiclassical approximation also the degeneracies 
are preserved and the levels of the semiclassical approximation have the same degeneracy as the exact 
ones. Therefore it is meaningless to compare the semiclassical error to the mean level spacing that is 
inversely proportional to the total number of energy levels. It is much more meaningful to consider 
the error within each group characterized by some symmetry and therefore by an additional quantum 
number. Since the existence of quantum numbers in addition to energy and of caustics are generic for 
integrable systems we expect some of the results found in this work to be generic for such systems. 
It is expected that the semiclassical approximation will fail for well defined groups of states and the 
fraction of such states, rather than the mean semiclassical error is a good measure of the quality of 
the semiclassical approximation. To our understanding in general the density of badly approximated 
states may depend on energy, contrary to the situation in the present work. The characterization of 
this dependence should be the subject of further studies. 

The results of this work may be relevant for mixed systems, for quantum states that are localized 
in the regular regions. For chaotic systems and for the chaotic component of mixed systems, there 
are no caustics and the wave functions usually spreads over all the chaotic region. The energy is the 
only identity of an eigenstate. Therefore the results of this work may not be of direct relevance for 
such systems and the mechanism for the destruction of the semiclassical approximation is different 
there. Some corrections to the trace formula for chaotic systems may be of a similar origin. 



Alonso and Gaspard computed the correction to the Gutzwiller trace formulae for billiards [20 
The correction is complicated but some of its parts have geometrical meaning. One term is pro- 
portional to the sum J2i tt where Cj is the chord length. Another contribution is proportional to 
Y,i Rc l s 3 -i where R is the radius of curvature of the wall and is the angle between the orbit and 
the normal to the wall at the incident point . When a caustic is close to a hard wall the typical 
chord length will be small and also cos0. Thus the correction to the leading order contribution of 
periodic orbits will diverge and the semiclassical quantization may fail. Therefore we expect that 
semiclassical methods will give poor results for contributions of orbits that are adjacent to a hard 
wall. For chaotic systems, because of ergodicity a typical orbit will not have a large fraction of chords 
near the wall of a billiard therefore this correction will not be dominate. One can expect that in 
such systems the semiclassical error will fluctuate around some average without extreme deviations. 
As is clear from the present work the situation will be different for integrable and mixed systems. 
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APPENDIX A: LEADING ORDER QUANTIZATION 



In this Appendix we use the "old quantum theory" in order to obtain the semiclassical quanti- 
zation condition (0). The method uses the classical momenta for quantization of action integrals. 
This result is not new but it may help to clarify the form of the semiclassical limit that was used 
in the systematic WKB expansion. Our system is a D dimensional spherical billiard and it is con- 
venient to describe it using a Hyperspherical coordinate system. The transformation between the 
Hyperspherical and Cartesian coordinates is: 

x D = r cos£ D _i 

= rsin£ D _iCos£ D _ 2 
xd~2 = r sin sin £c_ 2 cos £x>_ 3 (Al) 



x 1 = r sin £ D _ X sin £ D _ 2 • • • sin £ 2 sin £1 
It is easy to prove that the Hamiltonian of a free particle in the spherical coordinates is given by 



+ 



+ 



PL-3 



r 2 sin ^ D _! r 2 sin £ D „x sin £ D _ 2 



+ ••• + 



Pi 



(A2) 



This Hamiltonian satisfies the Staeckel conditions [jnj], therefore Hamilton's characteristic function 
is completely separable 



W(q)=J2W l (q i 



The Hamilton- Jacobi equation in these coordinates is 



(A3) 



' 8Wr 

dr 



1 



+ 



'dW, 



+ 



D-l 



sin £ D _ 2 



'dW, 



sin 
1 (dWt 



'dW, 



\ 2" 



sin 2 6 V 



2mE. 



(A4) 



The brackets in this equation must be constants of motion since they depend upon different variables. 
Since the constants are positive we will define af as the value of the zth brackets. We can write the 
Hamiltonian as 



n 



+ 



2m 2mr 2 



(A5) 



where L = a^-i is the generalized angular momentum. We will calculate «d-i by calculating the 
actions related to the angular variables. The first coordinate is cyclic thus 



«Si = 27r«i 



All the other actions are similar and are computed in 
Si = i -^dii = 



\ 



aj r^^rd^i = 2ii(a>i - a^i) 



sin & 



(A6) 



(A7) 
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The WKB quantization is done by quantization of the actions Si = 27rh(li + 7«/4) where 7, is the 
Maslov index of the action. Now the angular momentum is given by 



h 1 



D 



L = <xd-i = i- = h (j^ih + h + h 

i \ i=2 / 

This is just the semiclassical angular momentum L sc . 
One additional quantization is needed in order to obtain the energy levels 



(A8) 



s r = <j> p r dr = 2 J r~~r y — 2^ r 



(A9) 



where R is the radius of the Hyper spherical billiard. Here the Maslov index is 3, since the contribution 
from the turning point is 1 and the hard wall contributes 2. A straight forward calculation leads to 
the quantization condition: 



Or 



V2mER 2 -L 2 -L 



L 2 



arccos 



2mER 2 



2nh{n + -). 



Introducing v and z, and identifying L = L sc , the equation takes the form, 



v 



\/ z 2 — 1 — arccos ( - 

z 



IT ( 1 

— — = 7T I // : - 



(A10) 



(All) 



that is just ( P2| 



APPENDIX B: CALCULATION OF THE SECOND ORDER CONTRIBUTION TO THE 

QUANTIZATION 



The objective of this Appendix is to calculate the contour integral fl24|). The contour must circle 
the turning point at r m j n and also the point r = R that represents the hard wall. It is convenient to 
compute this integral by taking the contour infmitesimally close to the real line. The parts of the 
contour that are parallel to the real r axis give the same real contribution, which can be computed 
directly. Calculating this integral will reveal that the classical turning point gives an infinite con- 
tribution. The part of the contour that circles the turning point also give infinite contribution, and 
the sum of all the terms is finite. 

The standard way to deal with integrals of this kind is to integrate and differentiate using a suitable 
parameter as a variable (usually the energy). Here one first integrates with respect to the energy. 
The number of integrations is such that the contribution from the turning point to the contour 
integral in the complex r plane converges. Then the contour integral is computed. The last step 
is to differentiate with respect to the energy to obtain the desired result. For this purpose ([23|) is 
written in the form 

= 3L 2 SC d 1 5Lj c d 2 1 1 fm , 

2 4mr 3 dE (2mEr 2 - L 2 sc y/ 2 24m 2 r 5 dE 2 (2mEr 2 - L 2 sc y/ 2 Sr^2mEr 2 - L 2 SC 

Now each term will be treated separately. Integration of the first one results in 

(B2) 



- 2mr3./,. m&2 _r. L 4mJPV - 2L , c l V 2mEIP 2L „ 2 
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Differentiation with respect to E and change of variable from E to z gives: 



§ = d b^-i + arccos ( v = db (t?=t + arccos <3 1 1 ' (B3) 



The next term can be treated in a similar manner 



R 5L^ C dr 



12m 2 r 5 ^2mEr 2 -L 2 c ' 
Substitution of y — - helps to perform the integral and to obtain 



(B4) 



5L 2 , — — — / 3 m E 1 

'2mER 2 - L 2 I J- — 



12m 2 v ac V4L2J2 2 4i? 4 



\ 5£ 2 


u 




arctan 


/ 8L SC 


I 2 



1 



2mER 2 7Y 



Lit,- 



(B5) 



This term has to be differentiated twice with respect to E, and then the variable E will is replaced 
by z. After some manipulations the contribution of this term is found to be 

8 2 h 5 1 5 1 5 /1\ 

+ — — arccos I - ] . (Bo) 



dE 2 4L SC UL SC (z 2 - l) 3 / 2 AL sr ^ : 

The third term is just 

f R dr 1/1" 

h = / nr- , = = 7^— arccos - ) . (B7) 

A/sffe Arj2mEr 2 -L 2 4L SC \zj 



And the final result is 

dh d 2 I 2 I I 5 

a * dr= dE- ^- J3 = 4Z^7P = T + 12L SC (, 2 -1)3/ 2 (B8) 

that reduces to f|24f). 



APPENDIX C: EXACT GENERALIZED ANGULAR MOMENTUM EIGENVALUES IN WKB 

The WKB series for the angular momentum in D = 3 dimensions were constructed by Salasnich and 
Sattin [T^|. Summation of this series leads to the correct form of the quantum angular eigenvalues, 
namely, (1 + 1)1. Here this result will be generalized to the angular momentum in arbitrary dimension 
D. 

The first step is to compute the form of the Laplacian in Hyper spherical coordinates. Then 
by separation of variables an equation relating the generalized angular momentum operator to its 
projection on a D — 1 dimensional space is found. We will assume that the Laplacian has the form 

+ (CD 

where A^, is the angular part that does not depend on the radial variable. We will use induction in 
the number of degrees of freedom to prove that the laplacian indeed has this form. First one should 
note that for D = 2 and 3 this assumption holds. Assume it holds for a D — 1 dimensional system 
and add an additional Cartesian coordinate xp. The Laplacian is now 
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A 



d 2 
d 2 %D 



A 



D-l 



Ao_i is the Laplacian in the space of the coordinates x\, x 2 , ■ ■ ■ , %d-i- Transforming x%, x 2 , 
to Hyper spherical coordinates and using the induction assumption we obtain 



d 2 d 2 D — 2 d 

d 2 xo d 2 r r dr 



— A s 



(C2) 



(C3) 



The transformation to Hyperspherical coordinates involves only the two coordinates r and xe>. All 
of the angular coordinates will not change and therefore neither will A| ) _ 1 . One substitutes 



x D 
r 



Rcos9 
Rsm9 



The computation of the laplacian is straight forward and leads to 

d 2 D 



A 



d 2 R 



l_d_ J_ 

R dR + R 2 D 



where 



d 2 



n 9 1 

2 cot 9— + — 
9^ sm^ 9 



2 /i^D-l 



(C4) 



(C5) 



(C6) 



is the D dimensional angular momentum operator. Since the systematic high order WKB expansion 
works only for ordinary differential equations we will assume that the exact form of the eigenvalues 
of the angular momentum in lower dimensions is known. The equation connecting eigenvalues of 
angular momentum of D and D — 1 dimensions is obtained by separation of variables 



T"{9) + (D - 2) cot(0)T'(0) 



(m + D — 3)m 
sir?(0) 



T(9) = -\ 2 T{9). 



(C7) 



where T(#) is the eigenfunction for the eigenvalue A 2 , while m(m + D — 3) is the eigenvalue of A| ) _ 1 . 
It is justified by induction in what follows (see (|U11|) ). This is the generalization of equation (4) 
The computation follows [18| with only minor changes. Instead of (5) in [18] one substitutes 



in 



T{9) 



F(0) 
(sin oy 



(CI 



where a 



D-2 



, . After some manipulations and substitution of x 

U 



+ | one obtains 



— F"( 



x) 



cos 2 (x) 



F(x) = EF(x), 



(C9) 



that is (7) of p~8[j . However, the values of the parameters are different. Here E = X 2 + a 2 and 
U = m(m + D — 3) — a(l — a). The WKB expansion of (|C9| ) is developed in [fT8| , p~3[] for any order. 
The quantization condition is given by a series which is then summed. This results in 



AU 2 



1 



(CIO) 



But here y/1 + AU 2 



2(m + ^2^)- Substitution of U and E leads to 
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A 2 = (n + m)(n + m + D - 2) (Cll) 

This is indeed the correct result for the generalized angular momentum eigenvalues. The semiclassical 
eigenvalue is found from ( |C10|) assuming E ~ A 2 , leading toA~n + m+ ^y^. Identifying I = n + m 
leads to (|X8|) . 
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